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ON THE INTERIOR CONSTITUTION OF THE EARTH AS RESPECTS 

DENSITY. 

By Dr. G. W. Hill, Washington, D. C. 

Nearly all the matter accessible to us is found to be porous. Thus the appli- 
cation of pressure to it tends to reduce the amount of porosity and, in conse- 
quence, augments the density of the mass. Moreover, the greater the pressure 
the greater is the increment of density. A familiar instance of this is the case of 
atmospheric air or a gas in which, provided the temperature remains constant, the 
density varies directly as the pressure. 

It is natural to think that the matter of which the earth is composed is not 
excepted from this law. At small depths, it is true, the rigidity of the earth's mass 
interferes with its exerting any pressure, as the existence of caves shows. But at 
great depths where the weight of the superincumbent mass becomes very great, 
it is extremely probable the molecular force of cohesion gives way in a manner 
which allows pressure to act ; which is illustrated by the behavior of ice in a 
glacier. 

I propose to see what conclusions we are led to by adopting this relation be- 
tween the density p and the pressure /, 

p = A-\-Bp. 

A and B are constants, A denoting the density at the surface, and B the rate of 
increase of the density per unit of pressure. In applying this formula to the at- 
mosphere and gases, we have by Boyle's law A = o. Let V denote the potential 
of the gravitating force of the whole mass, and let us neglect the effect of the cen- 
trifugal force arising from the rotation of the earth. Then pressure being sup- 
posed to act as though the whole mass were fluid, hydrostatics furnishes us with 
the equation 

dp = pdV. 

F being restricted to points on the surface or in the interior of the mass, it satisfies 
the partial differential equation 

The three equations now written may be regarded as determining the three 
unknowns p, p, and V. 

By the elimination of Vandp we get 
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It will be seen that the constant A has disappeared from this equation. By Boyle's 
law in the case of gases A=o; that is, the matter is capable of attenuating itself 
to an infinite degree, a thing very improbable. But the introduction of the con- 
stant term A, and consequent supposition of a limit to the attenuation, does not 
change the differential equation which p satisfies. This partial differential equa- 
tion contains the whole theory of gases under a uniform temperature contained in 
vessels of any figure, and acted on by any gravitating forces ; also the theory of 
atmospheres surrounding solid nuclei of density as heterogeneous as we please, 
and of any figure. The truth of the equation is not at all invalidated by any dis- 
continuity m p or B; these quantities may change the law of their values as often 
as the problem demands. 

The very simple integral of this equation in the case of the earth's atmos- 
phere, when the attraction of the atmosphere on itself is neglected, is well known. 
It is our object here to examine the special solutions of this equation which are 
defined by the equation, 

p = function [i/(^^ +y + ^)'\- 

In this case, making r = ^{x^ + y + ■s^), the partiid differential equation is re- 
duced to an ordinary one and becomes 

^ + ^Br'p = o. 



or, as it may be written, 



d^ (r log p) , „ 



To simplify this, let us put 

s = 47iBr^p. 

Then .s being made the dependent variable, we have 

, „ d . log s 

d . r^ 3-2- 

dr 

5 Y s — 2 = 0. 

dr 

And if log r = t^, it becomes 

(/^ log s , d . log .f , 



di? dv 

Furthermore, if — " — = u, this differential equation of the first order between 
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u and s is obtained 

du 2 — {u -\- s) 

ds us 

This being integrated, and « obtained in terms of s, or s in terms of «, r is 
given by the equation 

or by the equation 

.A 



in which K is an arbitrary constant. And, if in the first of these values of r, ^Bf^p 
is substituted for s, the equation will be obtained which determines /> as a func- 
tion of r. 

The differential equation in u and j is a particular case of the general form 

Pdx + Qdy = o, 

where P and Q denote algebraical functions of x and y of the form 

Ax^ + Bxy ■\- Cf ■\- Dx + Ey -\- F. 

Mathematicians have been able to obtain the integral of this, in finite terms, only 
when the constants A, B, etc. satisfy certain equations of condition.* Unfortu- 
nately, the differential equation under consideration does not belong to any of 
these particular cases. Recourse must be had to series or other methods of ap- 
proximation for the determination of the relation between u and s. However, the 
differential equation itself will furnish the properties of the family of plane curves 
it defines. 

Thus u and .f denoting the rectangular co-ordinates of a point in a plane, the 
differential equation gives immediately the means of drawing the tangent to the 
curve which passes through this point. Excepting at the two singular points 
whose co-ordinates are « = o, .f = 2 and u = 2,s = o, for which the expression 
of the tangent takes the indeterminate form 

du o 

ds o' 

the system of curves do not intersect each other, since there is but one value of 

^ for given values of u and s. Since the differential equation is satisfied by the 

condition s = o, the axis of u is itself one of the system of curves, and no curve 
can cross it except at the point « = 2. If, in the differential equation, we substi- 

*See Xj\ouy\l\e, Journal de Mathimaliques, 2e Series, Tom. Ill, p. 417. 
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tute 2 -\- du for u, and ds for s, it is clear that only one curve passes through this 

point, and that its tangent here is given 

by the equation dujds == — ^. The 

axis of u, between the points « = 2 

and « = 00 , is an asymptote to the 

whole system of curves. The axis of 

s is intersected at right angles by the 

system of curves. Investigating what 

occurs at the point s = 2 on this axis, 

we substitute du for u and 2 -{- ds for 

s, and obtain for determining dujds at 

this point the following quadratic 



m 



■+4^i 



o. 




the roots of which are imaginary. 

Hence no curve passes through this 

point, and it is easy to see that the system of curves makes an infinite number of 

turns about it. 

The tangent to any curve, at its intersection with the straight line whose 
equation is « + .r =: 2, is parallel to the axis of s. When u and s are both very 
great, the tangent to the curve approximates to parallelism with the axis of s^ 
When s is very great and u small in comparison, the differential equation becomes 
approximately 



du 
ds 



i; 



or integrated, 



«* = 2 (Jq — s), 



if jj is the value of s when u= o. Hence the curves in the vicinity of the axis 
of s approximate to the parabola, in measure as we recede from the origin of co- 
ordinates. 

It is very easy to draw the curves connecting all the points possessing parallel 
tangents. For convenience let a denote the common value of dsjuu for these 
points; then the differential equation furnishes 

(« + a) (j + «) = a (a + 2). 

Thus these curves are equilateral hyperbolas having their asymptotes parallel to 
the axis of co-ordinates. 

Thus much in regard to the properties of the curves defined by the differ- 
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ential equation under consideration. But, for the special physical problem we 
have in view, there is no necessity to attend to the course of the curves through 
the whole plane. The density being supposed to increase with augmentation of 
pressure, B is necessarily positive, and r and p, from the nature of the problem, 
being the same ; s is likewise a positive quantity. There is then need only of con- 
sidering the curves on the positive side of the axis of u. Moreover, since 

d . log r ft dr ' 

and dp/dr is always negative when the force is directed towards the centre of the 
mass, there is no need of attending to the curves in the portion of the plane for 
which « > 2. 

Before proceeding to the special problem we have in hand, I propose to illus- 
trate the general theory by considering the density of the earth's atmosphere. It 
must be remembered that, in the usual manner of treating this question, the 
attraction of the atmosphere on itself is neglected; here, however, it is taken into 
account. Boyle's law being supposed to hold exactly, we shall have 

p = B/>. 

To integrate the differential equation between u and s, it will be necessary to ob- 
tain from observation the initial values of these two variables which hold at the 
surface of the earth. Let us denote these by «„and s^; and by a similar notation 
the values of all the variables at the earth's surface. The values of % and s^ re- 
sult from those of certain well-known physical constants. 

Let Z> = the density of mercury, 

k = the altitude of the barometer, 
£■ = the force of gravity,. 
R =^ the mean density of the earth. 



From an equation just given we have 

-0 = ^0 [^^^] +2 



But we also evidently have 






p,=gDh, 



Kdr] 



gP<^- 
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Substituting these values, «o = 2 — ^j^ ■ 

Thus it is apparent that u is independent of the units assumed Jfor the measure- 
ment of lengths and densities. In the next place 

A ^^^' 

But we have £■ = ~ — - . ~^ = ^ ^Rr„ . 

3 ^0 

Thence we get j„ = A^KBr^p^ = ^)M. ' 

Thus .f is also independent of the just mentioned units. 

Let us adopt the following values of the constants which enter into the ex- 
pressions of «o and jpi — 

^0 = 6365419 metres, 
A = 0.76 metres, 
Po= 0.001293 187, 
I)= 13.596, 
R= 5.67. 

The value of jO„ is that found by Regnault* for the temperature 0° of the cen- 
tigrade scale and the given altitude of the barometer; r, is the distance of his 
observatory from the centre of the earth according to Bessel's dimensions of the 
terrestrial spheroid ; and the value of R is that determined by Baily in his repe- 
tition of the Cavendish experiment. With these data we obtain the following 
values of «„ and s^ : — 

»o== — 794-6425, 
s„= 0.5450835. 

Having these initial values we can easily integrate the differential equation 
connecting u and s by mechanical quadratures or series, in the direction of s 
diminishing until s becomes so small as to be of no account. The corresponding 
values of r and p could then be found as we have already explained. However, 
the differences between the numerical values obtained by this method and those 
resulting from neglecting the action of the atmosphere on itself would be in- 
sensible. 

We pass now to the problem of the mass of the earth. Let us here denote 
the values of the variables which hold at the centre by the subscript („). If the 

*Memoires de P Academie des Sciences de Paris, Tom. XXI. 
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density at the centre be finite we must have j, = o; and the differential equation 

ds us 

du 2 — {u -\- s) 

shows that u^ = 2, else s would be o for all values of u. Hence the curve we 
have to consider, in this case, is the single one which passes through the singular 
point u = 2,s = o. 

The mass included in the sphere whose radius is r, is 

r 



if 



_ I , d. log fi 
~ B dr 

= ^r{2—u). 

Hence, denoting the values of the variables at the earth's surface by the sub- 
script (1), and R denoting, as before, the mean density of the earth, we shall have 

4^p,.s n(2 — «i) 

3 ' ~ B ■ 

Whence we derive £ = ' ,, 7* ■ 

and -^1 = 3(2 — «,)§. 

Then if we draw in the plane the right line whose equation is 

^ = 3^(2 — «), 

the co-ordinates of its intersection with the curve defined by the differential equa- 
tion and passing through the singular point u = 2, s^=o, will be the values of u^ 
and s^. This right line passes through the point « = 2, j = o, and it is readily 
ascertained from the differential equation that upon this curve u constantly dimin- 
ishes as s augments until it becomes o. The lines can therefore intersect on the 
positive side of the axis of s only when 

where OC is the distance from the origin of the point where the mentioned 
curve crosses the axis of s. 

In order to illustrate the general theory by an application, I have computed 
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by mechanical quadratures the values of the variable s and the function necessary 
for obtaining r. For this purpose it will be well to substitute for the independent 
variable u the variable s ^ 2 — u. The results obtained are given in the follow- 
ing table at intervals of o.i in ^'i — 



z 


J 


siz 


/- dz 
J t — g 


log r 


log sjr" 


0.0 


0.000 


3.000 


00 


— 00 


0.4771 


O.I 


0.294 


2.940 


— 1.1360 


9.5065 


0-4553 


0.2 


0.576 


2.879 


- 0.7737 


9.6640 


0.4323 


0.3 


0.846 


2.818 


- 0.5546 


9.7592 


0.4088 


0.4 


1. 103 


2-757 


- 0.3938 


9.8290 


0.3845 


0.5 


1-348 


2.695 


0.2646 


9.885 I 


0-3594 


0.6 


1.580 


2-633 


— O.I55I 


9.9326 


0.3333 


0.7 


1.799 


2.570 


— 0.0589 


9-9744 


0.3061 


0.8 


2.005 


2-507 


-f 0.0279 


O.OI2I 


0.2780 


0.9 


2.198 


2.442 


+ 0.1078 


0.0468 


0.2485 


I.O 


2.378 


2.378 


+ 0.1825 


0.0792 


0.2176 


I.I 


2-543 


2.312 


+ 0.2533 


O.I 100 


0.1854 


1.2 


2.695 


2.246 


+ 0.3213 


0.1396 


O.I5I4 


1-3 


2.832 


2.178 


+ 0.3874 


0.1682 


0.1 1 js 


1.4 


2-953 


2.1 10 


+ 04522 


0.1964 


0.0776 


1-5 


3.060 


2.040 


+ 0.5163 


0.2242 


0.0372 


1.6 


3-149 


1.968 


+ 0.5806 


0.2522 


9-9939 


1-7 


3.222 


1.895 


+ 0.6457 


0.2804 


9-9473 


1.8 


3-276 


1.820 


+ 0.7123 


0.3094 


9.8966 


1.9 


3-310 


1.742 


+ 0.7816 


0.3394 


9.8414 


2.0 


3-322 


1. 661 


+ 0.8547 


0.3712 


9.7791 


2.1 


3-309 


1-576 


+ 0.9336 


04055 


9.7088 


2.2 


3.265 


1.484 


+ I.O2I5 


0.4436 


9.6266 


2.3 


3-182 


1384 


+ 1. 1239 


0.4881 


9-5365 
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Let us suppose that the surface density of the earth p^ = 2.7 and the mean 
density R = 5.67. Then at the surface of the earth the value of slz must be 



Pi. 
' R 



1 .4286. 



By interpolating in the table it is found that this value corresponds to the follow- 
ing values of tne principal variables: — 

2=2.257, 
s = 3.224, 
log r = 0.4681, 

log p = 9-5722. 

Now the last two quantities are the logarithms of the surface values of the radius 
and the density measured in such units as in every case will give the simplest 
values to the arbitrary constants. But let us take the radius at the surface as the 
linear unit, and represent the surface density as 2.7. Then to reduce the num- 
bers so as to correspond to these units, it is evident we must add 9.5319 to the 
logarithms in the column of log r, and 0.8592 to the logarithms in the column 
of log sjr^. Thus are obtained the. following corresponding values of r and p: — 



r 


P 


r 


9 


0.000 


21.69 


0.469 


10.25 


0.109 


20.63 


0.501 


9-43 


0.157 


19-57 


0.535 


8.65 


0.195 


18.54 


0.570 


7.88 


0.230 


17-53 


0.608 


7-13 


0.261 


16.54 


0.649 


6.40 


0.291 


15.58 


0.694 


5.70 


0.321 


14.63 


0.743 


5.02 


0.350 


13-72 


0.800 


4-35 


0.379 


12.81 


0.866 


370 


0.408 


"-93 


0.945 


3.06 


0.438 


11.08 


1. 000 


2.70 



It will be noticed that the density at the centre is almost double of that given 
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by Laplace's formula ; and it seems that this supposition as to the law of density 
will not fit the phenomena as well as the latter. 

The limit beneath which the ratio p^jR cannot be reduced without the prob- 
lem failing to have a solution, is of interest. If the curve employed for the solu- 
tion of this problem is prolonged until its tangent passes, through the singular 
point on the axis of u, which it plainly must do before the curve crosses the axis 
oi s a. second time, this tangent affords the limit sought for the ratio 3 PiJR. The 
tangents of the curves, at the points of the plane whose co-ordinates satisfy the 
equation 

2 — {u -\- s) U — 2 

US S ' 

pass through the mentioned singular point. This equation in a simpler form is 

.J=:(l -\-u){2—u), 

which consequently represents a parabola passing through both singular points 
and having its axis parallel to that of s. By the employment of mechanical quad- 
ratures, the following additional points of the curve have been obtained : — 



s 


z 


s 


z 


30 


2.420 


2.3 


2.499 


2.9 


2.458 


2.2 


2.478 


2.8 


2.486 


2.1 


2.446 


2.7 


2.505 


2.0 


2.403 


2.6 


2.515 


1-9 


2.345 


2.5 


2.518 


1.8 


2.264 


2.4 


2.513 


1-75 


2.204 



From these it is evident the point u= — 0.2, s= 1.76 which lies -on the just- 
mentioned parabola is dso very nearly on the employed curve. Hence if PilR is 
less than a fraction which is approximately ^, there is no solution. 

The number of solutions in any particular Ccise is deserving of notice. The 
integral 



f 



s — s 



is proportional to the value of log r. It does not become infinite until the curve 
has made an infinite number of turns about the singular point on the axis of J- 
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This may be shown by a transformation of variables. Let us adopt polar co- 
ordinates, the singular point being the pole, and thus put 

.f = w cos -\- 2, 
^ = w sin d -\- 2. 
The differential equation then becomes 

flfw w sin 6 cos* + sin* ^ -|- sih ^ cos d 



w w cos d sin* ^ + I + sin* d — sin ^ cos d 

And we have 

dz r dd 



dd. 



J s — 2 J ^ 



w COS d sin* ^ + I + sin* — sin ^ cos d ' 

The denominator of these expressions cannot vanish unless w exceed 2, and it is 
plain that it remains positive and finite for all values of 0. Thus r becomes in- 
finite only when d does. Consequently there are an infinite number of solutions 
when p^R = ^; and a less number when pJR is either less or greater than this. 
With the value we have attributed to this fraction in the case of the earth, the 
course of the curve shows that there is but one solution. 
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Derive geometrically the usual expressions for the radius of the circle in- 
scribed in a triangle, and for the area of the triangle, in terms of the sides. 

[JF. H. Loud:\ 
165 

In any triangle ABC let a circle be inscribed touching the sides AB, BC, CA 
in N, L, M respectively. Let the centre of this circle be joined to the vertices, 
and from let OP, OQ be drawn perpendicular respectively to OC, OB, and cut- 
ting BC'm P and Q. Then if NP and AQhe drawn, these lines will be parallel 
as will also AP and MQ. [F. H. Loud:\ 
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A CIRCLE cuts a parabola and the centroid of the four points of intersection 
is found. What is the locus of the centre of the circle if this point be fixed? 

{W.M. Thornton:\ 



